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ice is considered as a major crop due to its demand globally. Pakistan is famous throughout 
the world to produce export quality rice which have healthy contribution in boosting the 
regional economy. Rice plant requite plenty of water for its proper growth and development 
therefore, water conservation is significant to maintain water reserves for a sustainable future. 
The main objective of this study was to identify day-to-day availability of water in rice fields from 
Germination to Ripening (GTR) using Carnegie Ames Stanford Approach (CASA) model. CASA 
model incorporates real-time parameter e.g., temperature, pressure, extraterrestrial radiations, Leaf 
Area Index (LAD, vapor pressure and sunshine hours to compute net-shortwave radiations (Ras), net- 
longwave radiations (Rn), net-radiations (Rn), actual incoming radiations (Rso), sensible heat flux (H), 
ground heat flux (G.) and finally the water stress (W). The averaged values of Ra, Reo, Rns, Raiand H 
were computed as 206, 319, 178, 34 and 124 (wm) respectively for GIR. Total expected sunshine 
hours were 1584h but we could receive only 874 h during GTR due to “off and on” cloud activity. LAI 
and G, were observed in inverse relation to each other. Ra, H and Go were used to compute the 
availability of water during GTR and concluded that the experimental site had surplus of water in early 
growth stages of rice crop having “W>0.5”, while milky dough and ripening stages were observed 
under water deficit (W<0.5). The results computed through CASA model were in accordance with 
actual observations collected through field survey up to 86%. 
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INTRODUCTION 

Rice (Oriza stavia L) is accounted as widely grown crop after wheat, over vast lands 
throughout the world [1]. Irrigated rice is an essential part of food security|2]. Rice is known 
as water baby, which needs excess amount of fresh water throughout its life cycle for its proper 
growth and development [3]. Various human activities are responsible for fresh water turbidity 
for instance industrial discharge, solid waste and the domestic discharge [3]-[6]. Asian 
countries are contributing about 90% in global rice production with utilization of 91% of 
global fresh water. It is projected that 15 M Ha of irrigated land and 22 M Ha of dry rice area 
may face water scarcity by 2025 [7]. Fresh water is essential to achieve targeted yield therefore, 
water saving techniques are of great importance for a sustainable and green future. Many 
researchers have devised a number of techniques to manage water supply in optimized range 
to receive healthy yield [8], [9]. All these techniques are related to the sustainable use of water. 
The main objectives of researches, which were performed in context of crop-water 
relationship, were to increase net crop production with less amount of water. 

Numerous factors are responsible to enhance productivity by provision of water in 
reasonable quantity e.g., soil type, soil pH and electric conductivity of soil [10]. Sand is not 
considered suitable for rice crop because it is perfectly drained. However, sandy clay is well 
drained therefore, considered less suitable while clay is best for rice cultivation because it is 
imperfectly drained and hold water for long times [11]. The clayish land has highest water 
holding capacity, so we can save a large amount of water with changing soil type [12]. The 
suitable range of electric conductivity for rice fields is 0.75-1.50 and the pH of soil must be 
between 5.4 and 7.1 pH [10]. 

Computation of water stress is a component of CASA model that initially requires 
various solar fluxes affecting crop growth and development which are based on sunlight and 
termed as Light Use Efficiency (LUE). The LUE directly affect water intake capacity of the 
crops. CASA model has been proved efficient to map productivity at global scales by fixing 
solar radiations [13]. Every crop needs certain amount of light for proper growth and 
development [11]. LUE is a very basic property of a plant [14] which express heterogeneity. 
The light used by plants is dependent upon a number of factors e.g. availability of chlorophyll 
content, leaf age, light intensity and crop’s developmental stage. Solar zenith angle, shape of 
crop canopy, leaf inclination angle and the leaf area are considered most influencing factors 
while computing the light used by plants. Dissimilar species are distinguished by computation 
of their ight used throughout their life cycle [15], [16]. There are many methods commonly 
used to compute light use efficiency including Eddy Covariance method [17], Quantum 
efficiency Renocking method and the productivity inversion method [18]. Photo Reflection 
Index (PRI) [19] is widely used to compute the variations in the light used by various crops 
while site-based measurements had been taken commonly to compute the LUE in the previous 
decades. 

This research aims at investigating rice fields suffering from water stress by 
incorporating various environmental factor essential for growth and development of rice crop 
e.g., temperature, pressure, sunshine hours and actual vapor pressure. It also aims at 
computation of variations in various fluxes throughout the rice growth period which effects 
the crop growth and development. These fluxes include net shortwave radiations (Ras), net 
longwave radiations (Rui), net radiations (Rn), possible incoming radiations (Rx), sensible heat 
flux (H), ground heat flux (G.) and actual incoming radiations (R.). CASA model incorporates 
Ras, Rat, Ra, Roo, H, Go and R, to provide water availability in rice fields on daily basis. 

Study site 

Hafizabad district was selected to execute this research. This district produces the rice 
of export quality. The study site is located at latitude between (31.4-32.2)°N and longitude 
between (73-73.5)°E. Monsoon directly hit the study site therefore, it receives large amount of 
water as rainfall which is good for rice cultivation. Water is distributed through connected 
network of water channels in an organized way. Weather is very severe where temperature 
vary from 0°C in winter to 48°C in summer. The study site is mapped in Figure 1. 
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Figure 1. Study site. 
The methodology to accomplish this research is described in Figure 2. 
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Figure 2. Flow of methodology. 

The basic environmental factors which were used to compute water stress in rice fields 
included daily fluctuation in temperature, possible sunshine hours, pressure, extra-terrestrial 
radiations, LAI, actual sunshine hours and the albedo of rice crop. 

The growth and development of crops is derived by sunlight [20] which is computed 
by estimation of Net radiations (R.). Ra shows the integrated impact of incoming shortwave 
radiations (R.;) and the earth’s emitted longwave radiations (Ru). Net Radiations effects the 
process of evapotranspiration therefore; it is important to compute R, to review agro-climatic 
interactions [21]. Mathematically, R, cab be computed as below [22], 

Rn = Rns + Roi (1) 


Where R,; ate net shortwave or incoming radiations, while Rai ate net longwave or 
outgoing radiations. Shortwave radiations are normally computed by incorporating albedo and 
actual incoming radiations [22]. 

Rns = (1 — a) Rs 


(2) 
Where « is the albedo of rice crop which is normally taken between a range of 0.15 - 
0.25 [23]—[25]and R, represents the possible incoming solar radiations. 
R, can be computed using the expressions as follows [27], 


Rs 


T T 
Rn = o[ mad mink (0.33 + 0.01LAI — 0.15,/éz)(0.84— + 0.15) 


Rso 3) 

Where o is Stephan Boltzmann’s constant and it’s value is 4.09% 10°MJ/k*M. Tinax 
and T'nin are daily fluctuations in temperature from GTR, R. is the actual solar flux while e, is 
the actual vapor pressure. We obtained Timin, Tmax,and e, from local metrological station. Extra- 
terrestrial radiations are represented by R,. R, is the solar flux which is recorded at the 
outermost boundary of our atmosphere [26]. R. depends upon solar angle that changes from 
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January to December. R, from GTR, was estimated using the website as below, 
http://www.ener.scu.edu/~emaurer/tools/calc_solar_cgi.pl. 

The actual solar flux that approached the earth’s surface by its interaction with clouds 
and other atmospheric contents are known as actual incoming solar radiations which were 
computed using the expression as below [22]. 

Rs = [As +Bs(—=)]Ra 
‘ (4) 

In equation 4, there are two calibration constants, A; and B, and the value of these 
constants is fixed as 0.25 and 0.50 respectively, N and n represents the possible sunshine hours 
and actual sunshine hours respectively. We acquired n, N from local metrological stations. 

If there is sunny day, even then about 25% radiations are scattered by their interactions 
with atmospheric gases and 75% approaches the surface of the earth. R,. can be computed 
according to Angstrom’s expression as below [28] 


Ro ae (A; + Bs) Ra (5) 


Ground heat flux (G.) 

Ground heat flux (G.) describes the penetration of heat in to the soil. It defines the 
heat conduction capacity of soil through which it is travelling that depends upon the 
concentration of vatious heat conducting elements present in a particular soil. G, actually 
control all bio-chemical reactions occurring in soil. The value of G, is very small therefore, it 
is ignored but it may create considerable errors sometime. G. was estimated using the 
expression as below [29], 


Go = Rn x (0.05 + 0.18 x exp (—0.521 x LAND) é 

Where LAI is the leaf area index and R, is the net radiations. Go is largely affected by 
LAI because an increase in LAI reduces the heat to reach the surface of earth and minimum 
leaf area at the time of germination leads to direct approach of sunlight in to the soil where it 
is grown. G, is minimum where the crop has maximum leaf area. In the study site, G. was 
observed minimum in the month of September where the leaf area reached at 6.65 m*m™. We 
obtained 20 field observations to record leaf area of the rice crop from July to October with a 
5-days temporal window as shown in Figure 3. 
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Figure 3. Temporal changes in LAI throughout the tice growth period. 

We fixed a 5 Days window between successive readings of leaf area because a leaf 
takes about 5 days to grow in considerable limits. Same temporal window was fixed by [11] to 
estimate the net rice production using CASA model. 

Sensible heat flux (A) 

Sensible heat is defined as the transfer of heat energy required for a plant to balance 
its temperature in accordance to its environment. It is difficult to compute the sensible heat 
flux accurately [30] therefore, it is taken as residual of energy balance [31].Bowen ratio (8) plays 
a vital role for computing H. Bowen ratio can be computed by taking a ratio of variations in 
vapor pressure and temperature of rice crop canopy using the expression as below, 

AT 
B = ¥z- 
Ae (7) 
Where y is a psychometric constant and its value can be computed using the expression as 
below [32] 
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y = 0.000665 x P 8) 


Variations in actual vapor pressure, temperature and pressure were acquired from local 

weather station and the sensible heat was computed using Bowen Ratio [31] as below, 
B 
H = (FP 5)(Rn~G) 3 

Where 8, Ra and G, are Bowen ratio, net radiations and the ground heat flux 
respectively. 

Water stress (W) determine the amount of available moisture content in soil. Water 
stress can be computed by the following expressions as below [33]. 
_R,-G,-H (10) 

R, = G, 

W vary between 0 and 1, where W=0 represents the oven dry soil and w = 1 shows the surplus 
of water [34] — [38]. 
Results and Discussion 
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Figure 4. Fluctuations in various solar fluxes throughout the rice crop life cycle. 

Figure 4 is showing fluctuations in different fluxes faced by rice canopy throughout 
the rice crop lifecycle. Rice germination started on June 20, 2019 and harvesting of the rice 
crop ended on October 20, 2019. During this period, the extraterrestrial radiations were 
declined due to change in solar angle. About 55640 wm” were received on the outer boundary 
of our atmosphere during 120 days of rice crop growth but we could actually receive a flux of 
27470 wm” that approached the rice crop canopy. The averaged estimates of Ra, Rso, Ras, Rat 
and R, were computed as 426 wm”, 319 wm”, 178 wm”, 34 wm” and 206 wm* throughout 
the rice crop life cycle. Figure 4 is showing many peaks and dips in Ra, Ras, and Rai that resemble 
with the shape of n/N mapped in Figure 5. The fluctuations in Figure 5 show that as the value 
of n/N approaches to 1, it represents a sunny day and its value near to zero represents a cloudy 
day. 
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Figure 5. Variations in actual and possible sunshine hours. 
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Figure 5 is showing that cloud activity throughout the life cycle of rice crop. We 
observed that the value of n/N must be near to zero in early growth stages of rice crop which 
express sever cloud activity and probably the rainfall, while it must be near to 1 in milky dough 
and ripening stages that results in sunshine situation which is good for ripening of rice crop 
properly. Figure 5 is showing that a total possible sunshine hour, which were expected to 
approach the rice crop canopy were 1584 h but we could receive only 874 h throughout the 
tice crop cycle due to cloudy situation on the study site. 
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Figure 6. Ground heat flux throughout the rice crop life cycle. 

Figure 6 was mapped by putting R, and LAI using equation 6 which is showing a 
decreasing trend. The gradual decline to G, curve is due to increase in LAI. As leaf area 
increases, it stops heat to penetrate and to approaches the ground. We computed fluctuations 
in ground heat flux and the results are shown in Figure 6. 

We applied a regression of second order to LAI and G, and found R*=0.8843 as shown 
in Figure 7. 
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Figure 7. A 2” order correlation between LAI and G.. 
We computed the sensible heat flux of rice crop by putting the values of Bowen ratio 
(6), Ra and G, in Equation (9) and the results are shown in Figure 8 as below 
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Figure 8. Fluctuations in H throughout the rice crop life cycle. 
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The shape of curve in Figure 8 resembles with the shape of R, curve in Figure 4. It is 
due to the reason that H and R, are in direct relationship, as R. increases, it need more heat 
for a plant to maintain its body temperature therefore, H increases. Sensible heat flux 
throughout the RLC was recorded as 124 wm” as averaged. We computed water stress in rice 
fields on daily basis using G., H and R, in Equation 10 and the results are shown in Figure 9 
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Figure 9. Water stress throughout the rice crop life cycle. 

Figure 9 shows that day where curve descends to zero below 0.5, it represents 
deficiency of water in rice fields, in comparison to the day where the curve crosses 0.5 toward 
1. The study site had excess of water (W>0.5) in germination, leaf emergence, tillering and 
panicle primordia initiations while the literal stages of growth were in water stress (W<0.5). 
On visiting the fields, we found that ripening stage of rice crop is water free stage while initial 
stages need penalty of water for proper development. Finally, we observed that all crop growth 
factors Ra, Ras, Ra, Go, H and W ate largely affected by the cloud activity (n/N). The curve in 
figure 9 is representing, if w<0.5, there is water stress and vice versa. 

Conclusion: 

Both monitoring and evaluations of environmental parameters (1.e., temperature, 
pressure, actual vapor presser and extraterrestrial radiations) play a vital role to understand the 
sequential growth stages of rice crop for optimum water use. One can determine the deficiency 
or penalty of water in rice fields by incorporating real-time ecological parameters. CASA based 
methodology was adopted to accomplish this research which is handy to apply at global scales 
to monitor the rice crop in a better way. 
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